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(71) We, National Research Devel- 
opment Corporation, a British corporation 
established by Statute, of Kingsgate House, 
66/74 Victoria Street, London, S.W.I., do 

5 hereby declare the invention, for which we 
pray that a patent may be granted to us, and 
the method by which it is to be performed, 
to be particularly described in and by the 
following statement: — 

10 This invention relates to a quadrature 
spatial frequency Fourier analyser. 

Fourier analysers find application, for ex- 
ample, in the fields of radiography and the 
testing of lenses and other optically trans- 

15 missive elements or systems, where it is de- 
sired to determine the optical transfer func- 
tion and/or line spread function of a light- 
transmitting element or system. By using an 
optical Fourier transformation a non-periodic 

20 and discontinuous function may be repre- 
sented in terms of a continuous function, e.g. 
as a graphically presented curve. For example, 
in the field of radiography one can measure 
the optical transfer function of an X-ray 

25 fluorescent screen photoelectrically by an area 
type masking method which can be applied 
easily to the fluorescent screen. 

The present invention concerns an analyser 
which, for a sequence of discrete spatial fre- 

30 quencies, can determine the amplitude of the 
one-dimensional Fourier transform of a finite 
area of an image, i.e. the image being regarded 
as a spatial distribution of light intensity. 
Known Fourier analysers can be divided 

35 into two distinct categories. Firstly, analysers 
dependent on scanning methods, wherein the 
optical intensity distribution along the scanned 
path is converted into an electrical analogue 
which is then analysed by suitable electronic 

40 devices to give a spectrum representative of 
power per unit bandwidth versus frequency. 
From a knowledge of the scanning speed this 



can be equivalently expressed as relative power 
per unit bandwidth versus spatial frequency. 

Such scanning methods present almost in- 45 
superable problems when trying to analyse an 
image which itself has been produced by a 
scanning technique. For example, noise 
measurement by photographic integration fol- 
lowed by scanned analysis suffers from limited 50 
accuracy, dynamic range and integrating time, 
due to photographic limitations, largely attri- 
butable to the non-linearity of the process. 

The second category of analysers are those 
using non-scanning methods, which have how- 55 
ever so far been restricted to analysing fine 
spread functions. In such methods the object 
is imaged on to a mask whose effective trans- 
mission varies sinusoidally in one direction. 
When a maximum of mask transmission coin- 60 
cides with the point of maximum intensity 
of the spread function, the light flux passing 
through the mask will be a maximum and will 
yield the cosine Fourier transform amplitude 
at the spatial frequency of the mask. If the 65 
mask is displaced through 90°, the light flux 
will give the sine transform amplitude. For a 
symmetrical spread function the sine transform 
vanishes. The phase angle relative to the maxi- 
mum intensity can be obtained as a function 70 
of spatial frequency from the ratio of the two 
transforms, and the curve representing the 
amplitude of the power versus spatial frequency 
can be found for each frequency from the 
quadratic sum of the two transforms. This 73 
method cannot however be applied to a 
dynamic image whose intensity distribution 
varies in time, because the phases of the com- 
ponents cannot be identified. 

The analyser of the present invention on 80 
the other hand differs from analysers of said 
second category in that it can measure phase 
relative to the mask position, which is assumed 
to be fixed relative to some fixed origin, where- - 
as such known analysers determine phase rela- 85 
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tive to the line spread maximum intensity. 
Hence the analyser of the present invention 
can determine the Modulus of the Fourier Co- 
efficient at a given spatial frequency regardless 

5 of the phase relationship between object and 
mask intensity distributions. 

In accordance with the present invention, 
there is provided a Fourier analyser for 
measuring the one-dimensional spatial fre- 

10 quency spectrum of an image comprising 
optical means producing multiple images of 
the intensity distribution of said image under 
analysis, filter means having transmission char- 
acteristics varying sinusoidally in one direction 

15 and arranged to provide simultaneously from 
said multiple images outputs indicative of the 
sine and cosine transforms of the image under 
analysis, and means arranged to quadratically 
sum said outputs and provide a unidirectional 

20 signal representative of the spectrum amplitude 
at the spatial frequency of said filter means. 

A feature of the present invention resides 
in the ability of the analyser to deal with 
images whose intensity distribution varies 

25 temporally as well as spatially. This means, in 
particular, that it can analyse noise distribu- 
tions and images generated by scanning 
methods. 

The analyser of the present invention has a 
30 number of unique facilities: it will give, for a 
range of discrete spatial frequencies, the ampli- 
tude of the Fourier transform of 

a) images whose Fourier components are of 
unknown phase; 
35 b) dynamic images, typically produced by 
a quantum-limited scintillation process; 
and, 

c) scanned images, which may or may not 
additionally belong in b). 
40 These potentialities result from the analyser 
including filters to give simultaneously the 
values of the sine and cosine transforms of an 
image and a means of simulating negative 
transmission, thus enabling the analyser to 
45 make band-pass rather than low-pass measure- 
ments, as well as generating a D.C. output 
signal, leading to simple processing methods 
and making quasi-infinite integration possible. 
In one embodiment of analyser according 
. 50 to the present invention, the object is simul- 

90 



taneously imaged on to two sinusoidal masks 
of relative phase displacement ^/2, the quad- 
ratic sum of the outputs then giving the spec- 
trum amplitude at the spatial frequency of 
the mask. The spatial frequency component 55 
being measured can be visualised as a unidi- 
mensional sinusoidal variation in intensity of 
the same spatial frequency as the mask, so 
that for a relative mask-image phase difference 
of zero the intensity and transmission maxima 60 
coincide. However, the value of the quadratic 
sum is independent of the relative phase angle, 
which consequently need not be known. 

In other words, the spatial frequency 
Fourier transform can be written 
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F(s)=/a(w) cos (ws) dw+/b(w)sin(ws)dw (1) 
= / S(w) cos (ws+£) dw (2) 

where 

F(s) is the spatial intensity distribution; 
W is the angular spatial frequency; 
a(w) is the cosine Fourier transform ampli- 
tude; 

b(w) is the sine Fourier transform ampli- 
tude; 

S(w) is the modulus of the complex angular 
spatial frequency spectrum; and 

^(w) is the relative phase angle at frequency 
W. 

The main output required from the analyser 
is a plot of S(w) versus W. 

From(l) and (2): 
S(w) 2 =a(w) 2 +b(w) 2 ; 



If the mask transmission maximum and 85 
spatial frequency component amplitude maxi- 
mum do not coincide but differ by an angle 3, 
this is equivalent to adding a constant error 
angle £ in equations (1) and (2), which become; 



F(s)=/a(w) cos (ws+S)dw+ /b(w)sin(ws+S)dw 
=/S(w) cos (ws+S+*)dw. 



(la) 
(2a) 



Equation (3) is unchanged under these circum- 
stances, as can be seen by inspection. 
Equation (4) becomes 



T -b(w) 1 
L a(w) J 



(4a) 



Thus the sole result is that the reference 
angle, relative to which p(w) is measured, has 
been changed by 8. Although the reference is 
quite arbitrary, and the only interest in phase 
concerns phase change with frequency and 
never absolute phase, the relative shift <&(w) 
cannot be measured unless the error S can be 
separated from the sum <5+f(w). For a dynamic 
image this will involve some form of instan- 
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taneous servomechanism to continually adjust 
the spatial position (phase angle) to maximise 
the output from the cosine transform filter. 
For a static image such separation could be 

5 attempted, the phase shift required to maxi- 
mise the cosine output being 8. 

The mask can have a band-pass property- 
only if both positive and negative transmis- 
sion is possible. This can be simulated for 

10 example by subtracting from the orthogonal 
transform signals a mean level signal, also 
derived from the object. 

The electrical output is unidirectional: 

1) For an unscanned static image the out- 
15 put is a steady current; 

2) for an unscanned dynamic image the 
output is a varying current. The total charge 
integrated electrically in a given time of ob- 
servation is the same as would result if the 

20 dynamic image cou ld b e optically integrated 
for the same time and the integral presented 
to the mask. Thus a difficult technical pro- 
blem in image storage is replaced by a simple 
electrical integration; 

25- 3) for a scanned image the output is a vary- 
ing current whose integral over one scan period 
gives the same charge as would be unscanned 
original image if displayed for one scan dura- 
tion. 

30 The form and positioning of the optical 
components and the form of the filters used 
may take several alternative forms. The accom- 
panying drawing illustrates one embodiment 
of analytical apparatus in accordance with the 

35 invention. 

As shown in the drawing, the analyser basic- 
ally comprises a mirror tunnel formed by two 
parallel mirrors Ml and M2. An original 
image 01 having an intensity distribution to 

40 be analysed is arranged to be produced at one 
end of the tunnel and an objective lens L is 
positioned adjacent to the other end of the 
tunnel. Virtual images VI of the intensity dis- 
tribution to be analysed are formed by the 

45 parallel mirrors Ml and M2 and three images 
are thereby formed in the plane of three spaced 
masks SI, S2 and S3, which may for example 
be formed on a single photographic plate. 
Behind each mask is a photomultiplier PI, P2, 

50 P3 or other linear photosensitive device pro- 
viding outputs A, B and C respectively. 

Masks SI and S3 are of sinusoidally vary- 
ing transmission in one direction and are dis- 
placed through jt/2 relative to each other. 

55 Mask S2 has a transmission characteristic 
which enables the mean luminance to be 
determined. 

The orthogonal transform amplitudes are 
obtained from difference circuits 10 and 11 

60 by subtracting the mean output B of photo- 
multiplier P2 from the outputs A and C of 
photomultipliers PI and P3 respectively. The 
orthoganol transform amplitudes are then 
quadratically added by suitable circuitry 12, 



for example, comprising squaring circuits and 65 
a summing amplifier. 

The optical duplication which is provided 
by the mirror tunnel Ml, M2 may alterna- 
tively be provided by using a prismatic beam 
splitter, or a large collimating lens and a 70 
number of smaller imaging lenses, for ex- 
ample. 

In an alternative embodiment four masks 
may be used, one pair containing alternate 
half cycles of a sine wave. Behind each mask 75 
of this pair a photomultiplier is again used to 
provide an output signal and the difference in 
electrical outputs gives the sine transform 
amplitude. A corresponding arrangement for 
the other pair of masks gives the cosine trans- 80 
form amplitude. 

In either of these embodiments the masks 
may be either of variable density type or of 
variable area type. 

Variable density masks vary sinusoidally in 85 
transmission in one direction and have con- 
stant transmission at right angles to it 

Variable area masks consist of an opaque 
background with a sinusoidal variation in trans- 
parent area in one direction. If a cylindrical 90 
optic is added to a variable area mask in the 
imaging path to defocus the image at right 
angles to the direction of mask variation, the 
mask will have the same effect on the light 
flux transmitted as in a variable density mask, 95 

In general, the analyser system can measure 
the one-dimensional spatial frequency spec- 
trum in any image, including those formed 
by a scanning process. 

' For a stationary image no temporal pro- 100 
blems exist but for an image which is quan- 
tum limited the integrated intensity distribu- 
tion will be dependent on the build-up time 
allowed. In such cases photographic integra- 
tion might seem the obvious method but the 105 
number of processing techniques involved, each 
capable of causing a considerable change in 
the recorded or reproduced image, and the 
non-linearity of the photographic process and 
its limited dynamic range makes this of 110 
limited use. As long as the quantised imaging 
process is truly stationary in time, so that 
sequential measurements each for the same 
statistically long time will yield the same 
results, the apparatus hereinbefore described 115 
will perform a transform analysis. In prin- 
ciple, by increased optical complexity, measure- 
ments throughout the spectrum could be made 
simultaneously, but this would be a unique 
measure of the spectrum during that time 120 
only and would have no general applicability 
unless the process were stationary in time, in 
which case the much simpler sequential method 
is just as valid. 

At any instant the image has a one dimen- 125 
sional Fourier transform whose amplitude 
varies with frequency. If the quadratic addi- 
tion of sine and cosine components is per- 
formed as a continuous process, a continuous 
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record of the variation with time of transform 
amplitude is obtained. However, in general, 
one seeks the true or "infinite time" spectrum, 
which can be obtained as the average of the 
vector indicated at any instant in the con- 
tinuous record or much more simply by integ- 
rating the sine and cosine signals for the 
required period and then quadratically adding. 

For truly random noise the average value 
of the sine and cosine transforms will be iden- 
tical, and so the phase difference obtained 
from a given ratio of integrated sine and 
cosine signals can be used as a measure of 
spatial randomness or lack of correlation in 
the image. The sensitivity of the test for 
randomness will depend on the inherent signal- 
to-noise ratio of the measuring apparatus. 

For the measurement of spectra using long 
integration times the technique is basically through 90° relative to each other and the 



outputs and provide a unidirectional signal 
representative of the spectrum amplitude at 
the spatial frequency of said filter means. 

2. A Fourier analyser as claimed in claim 
1, wherein said filter means act as band-pass 
filters. 

3. A Fourier analyser as claimed in claim 
1 or 2, wherein said filter means comprise 
variable density masks. 

4. A Fourier analyser as claimed in claim 1 
or 2, wherein said filter means comprise vari- 
able area masks. 

5. A Fourier analyser as claimed in any 
preceding claim, wherein said optical means 
provides three images of the original, wherein 
said filter means comprises three masks, two 
of which are of sinusoidally varying trans- 

' one direction and are displaced 



D.C. pulse integration and involves no high- 
frequency analogue circuits. For a continuous 
record, analogue computation at frequencies 
as high as 20Mc/s is possible. 

The spatial frequency range over which 
the apparatus will function is determined at 
low frequencies by the bandwidth of the filter- 
ing process (e.g. size of mask) and at high 
frequencies by the quality of optics obtainable, 
the geometrical and mechanical design chosen, 
and by the quality of the sinusoidal masks or 
grids. 



third of which has a transmission character- 
istic enabling the mean luminance to be de- 
termined, wherein linear photosensitive devices 
are positioned at the output side of each mask, 
and wherein the orthogonal transform ampli- 
tudes are obtained by subtracting the mean 
output of said third mask from the outputs of 
the respective other two masks. 

6. A Fourier analyser as claimed in any of 
claims 1 to 4, wherein said optical means 
provides four images of the original, and 
wherein said filter means comprises four masks, 



Within limits set by the physical design of the sine transform amplitude being obtained 



the apparatus the bandwidth of measurement 
can be shaped as required, the actual band- 
width being reciprocally related to the length 
of the sinusoidal mask or grid and the shape 
determined by the transform of the envelope 



angular envelope of a simple truncated sinu- 
soid. By suitable choice of optical reduction or 
magnification any required frequency range 
can be obtained. 

WHAT WE CLAIM IS: — 

1. A Fourier analyser for measuring the one- 
dimensional spatial frequency spectrum of an 
image comprising optical means producing 
multiple images of the intensity distribution 
of said image under analysis, filter means 
having transmission characteristics varying 
sinusoidally in one direction and arranged to 
provide simultaneously from said multiple 
images outputs indTcative of the sine and 
cosine transforms of the image under analysis, 
and means arranged to quadratically sum said 



from the difference in the outputs of two of 
said masks, and the cosine transform ampli- 
tude being obtained from the difference in 
the outputs of the other two masks. 

7. A Fourier analyser as claimed in any of 
claims 1 to 5, wherein said optical means 
comprises a mirror tunnel formed by two 
parallel mirrors having the image under 
analysis located therebetween. 

8. A Fourier analyser as claimed in any of 
claims 1 to 6, wherein said optical means 
comprises a prismatic beam splitter. 

9. A Fourier analyser as claimed in any 
preceding claim when used for analysing 100 
dynamic images, which includes a servo- 
mechanism arranged to continually maximise 

the output of the cosine transform filter means. 

10. A Fourier analyser substantially as here- 
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inbefore described with reference t< 
panying drawings. 

POLLAK, MERCER & TENCH, 
Chartered Patent Agents, 
Audrey House, Ely Place, 

London E.C.I. 
Agents for the Applicants 
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